The age incidence of any cancer can be explained by a one-mutation model 
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kjr)! Abstract. We propose a one mutation model for cancer with a mutation rate that 

increases with time. Under rather general hypotheses the number of mutations is neces- 
sarily a (non homogeneous) Poisson process with the prescribed mutation rate. We show 
that the cumulative probability of cancer up to time t is, up to a multiplicative constant, 
an antiderivative of the mutation rate. 
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Introduction Mathematical models explaining age incidence of cancers go back to at 
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least Armitage and Doll (1954). The statistical data show that age incidence follows power 
laws for several different cancers (with different powers for different cancers). This lead 
to the hypothesis that a given cancer arises after several successive mutations and that 
there is a simple relation between the number of necessary mutations and the power law of 
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the incidence. This hypothesis has yet to be proved biologically and several authors have 
argued that cancers could be modelled by a two stage process, see Armitage (1985) and 
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Moolgavkar and Knudson (1981). Recently, Michor et al. (2006) proposed a one mutation 
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model to explain the age incidence of chronic myeloid leukemia which increases as a third 



power of age. In this paper, we show that in fact any power can be explained by a one 
mutation model. 

In this note, we propose a one mutation model with a mutation rate that increases 
with time. The number of mutations is necessarily a (non homogeneous) Poisson process 
with the prescribed mutation rate. We show that the cumulative probability of cancer up 
to time t is, up to a multiplicative constant, an antiderivative of the mutation rate. 

The model We assume that a given cancer appears after the following two steps. 
First we must have a certain mutation and then a clonal expansion of the mutant. Let 
N(t) be the number mutations up to age t. We make the following hypotheses for N. 

(i) N(0) = 0. 
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(ii) The number of mutations in disjoint time intervals are independent random vari- 
ables. That is, for any integer k > 1 and any sequence of times t < ti < . . . < tk we have 
that the random variables N(t n ) — N(t n -i), 1 < n < k, are independent. 

(iii) The probability that there is exactly one mutation during the time interval (t, t+h) 
is of the order of /j,(t)h as h goes to 0. The mutation rate \x > is a function of time. More 
precisely, we have 

lim P Wt + ft )-iv( t ) = i) 

h^o h 

(iv) The probability that there are two or more mutations during the time interval 
(t, t + h) is of order less than h as h goes to 0. That is, 

Um PjNjt + h)- N(t) > 2) = Q 

h-»0 h 

These four hypotheses are quite natural for a mutation process. They model the 
complete randomness of mutations. Under these hypotheses the process N is necessarily 
a Poisson process with rate 

m{t) = I ji{s)ds. 



o 



See, for instance, Section 5.4.1 in Ross (2006). In particular, 

Tn(t) 

P(N(t) = k) = exp(— m(t)) — — — for any integer k > 0. 

K. 

In order to provoke a cancer a given mutation must be successful. This might mean, for 
instance, replacing the healthy cells of an organ by the mutated cells or a clonal expansion 
of the mutated cells. Let a(t) be the probability that a mutation appearing at time t will 
expand. Assume that each mutation expands or not independently of all other mutations 
and of the process N. Under this assumption the number of mutations that will expand, 
denoted by C(i), is also a Poisson process with rate 



M(t) = [ a(s)n(s)ds. 
Jo 



For this fact, see (5.9) in Kingman (1993). Let I(t) be the probability of cancer by time 
t. It seems reasonable to assume that the time it takes for a mutation to expand in the 
body is small (weeks or months) compared to the time it takes for the appearance of a 
successful mutation (years or decades). Hence, we may approximate I(t) by P(C(t) > 1). 

I(t) = P(C(t) > 1) = 1 - exp(-M(t)). 
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Even at old age the probability of a given cancer is small. Therefore, at all ages t, M(t) 
must be close to and we may approximate 

I(t) ~ M(t) = [ a(s)fi(s)ds. 
Jo 

Suppose now that a is a constant (the probability of success for a given mutation does 
not depend on age). Then the cancer incidence I is up to a multiplicative constant an 
antiderivative of the mutation rate fi. Hence, in order to model a given incidence I one 
may define fx as constant times the derivative of I and then adjust the constant. 
For instance, if I(t) = ct 1 then set 

cy 

fi(s) = /Uos 7_ where fio = — . 

a 

There are several possible interpretations for a. It may be the fixation probability 
(the probability that eventually all cells of the organ will have the mutation) of a mutation 
in a given organ. Using a Moran model as in Michor et al. (2006) the fixation probability 
is given by 

1 - 1/r 

& = ; 

1 _ l/ r n 

where r is the relative fitness of the mutation and n is the number of cells in the organ. 

We get another expression for a if we assume that a mutant expands according to a 
branching process. Assume that each mutant cell divides with probability p or dies with 
probability 1—p independently of each other and of everything else. Starting with a single 
mutant the probability of an ever lasting expansion is 

a = 2 - - for p > 1/2. 
p 

This expression for a is obtained by observing that this is a branching process and 1 — a 
(the probability that the mutation dies out) is a solution of 

x = (1 — p) + px 2 . 

See, for instance 1.9 in Schinazi (1999). 

Discussion. We point out that plausible hypotheses (i)-(iv) the number of mutations 
in a one-mutation model is a non-homogeneous Poisson process. By neglecting the time of 
fixation (or expansion) and the time of detection the age incidence for this model is (up to 
a multiplicative constant) an anti-derivative of the mutation rate. Note that we would get 
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the same result if instead of neglecting these two times we would take them to be constant 
(as proposed by Moolgavkar and Knudson (1981)). According to our result if the mutation 
rate is constant the incidence is necessarily linear in age. In contrast, Michor et al. (2006) 
claim that the incidence may increase as the cube of age for a constant mutation rate in 
a one-mutation model. They, however, do not neglect the time for a mutant to fixate nor 
the time to detect the tumor. They take the first time to be constant and the second one 
to be random. 

The idea of using a non-homogeneous Poisson process for carcinogenesis goes back to 
at least Kendall (1960) although it seems to have been ignored by the recent literature . See 
also Moolgavkar and Venzon (1979), Moolgavkar and Knudson (1981)) and Whittemore 
and Keller (1978). Our main result, however, seems to be new. 
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